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^ ■ Abstract 

In a frameworli that makes close contact with modern effective field theories for non-relativistic 

bound states at zero temperature, we study the real-time evolution of a static quark-antiquark 

Q ' pair in a medium of gluons and light quarks at finite temperature. For temperatures ranging from 

values larger to smaller than the inverse distance of the quark and antiquark, 1/r, and at short 

, distances, we derive the potential between the two static sources, and calculate their energy and 

I thermal decay width. Two mechanisms contribute to the thermal decay width: the imaginary part 

G\ ■ 

• of the gluon self energy induced by the Landau damping phenomenon, and the quark-antiquark 



X 



color singlet to color octet thermal break up. Parametrically, the first mechanism dominates for 



o 

00 

I temperatures such that the Debye mass is larger than the binding energy, while the latter, which 



we quantify here for the first time, dominates for temperatures such that the Debye mass is smaller 
than the binding energy. If the Debye mass is of the same order as 1/r, our results are in agreement 
with a recent calculation of the static Wilson loop at finite temperature. For temperatures smaller 
than 1/r, we find new contributions to the potential, both real and imaginary, which may be 
relevant to understand the onset of heavy quarkonium dissociation in a thermal medium. 

PACS numbers: 12.38.-t,12.38.Bx,12.38.Mh,12.39.Hg 
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I. INTRODUCTION 



The study of heavy quark-antiquark pairs in a thermal medium at temperature T has 
received a lot of attention since it was suggested that quarkonium dissociation due to color 
screening may be a striking signature of the quark-gluon plasma formation [1] . Based on this 
idea and the assumption that medium effects can be understood in terms of a temperature 
dependent potential, the problem of quarkonium dissociation has been addressed in terms 
of potential models with screened temperature-dependent potentials over the past 20 years 
(see e.g. Refs. 2, for some representative works). A derivation from QCD of the 

in-medium quarkonium potential has not appeared in the literature so far and expectedly 
not all medium effects can be incorporated into a potential. A first step towards a QCD 
derivation of the quarkonium potential at finite temperature has been a recent calculation 
[sl of the static Wilson loop in the imaginary-time formalism at order ctg. After analytical 
continuation to real time, the calculation shows a real part, which is a screened Coulomb 
potential, and an imaginary part that may be traced back to the scattering of particles in 
the medium carrying momenta of order T with space-like gluons, a phenomenon also known 
as Landau damping. Some applications can be found in jy, Q] • First principle calculations 
of quarkonium properties at finite temperature include calculations of Euclidean correlation 
functions in lattice QCD and the reconstruction of the corresponding spectral functions 
using the Maximum Entropy Method. At present, however, a reliable determination of the 
quarkonium spectral functions from the lattice data appears very difficult due to statistical 
errors and lattice discretization effects (see discussion in Ref. 8| and references therein). 

In this work, we will study static quark-antiquark pairs in a thermal bath in real-time 
formalism (see e.g. jlj]) and in a framework that makes close contact with effective field 
theories (EFTs) for non-relativistic bound states at T = 1^. In this framework, we will 
address the problem of defining and deriving the potential between the two static sources, 
and we will calculate their energy and thermal decay width. We will describe the system for 
temperatures that range from larger to smaller values with respect to the inverse distance 
of the quark and antiquark, 1/r. In some range, we will agree with previous findings; 
for temperatures lower than the inverse distance of the quark and antiquark, we will find 
new contributions to the potential, both real and imaginary, with a non trivial analytical 
structure. In particular, we will point out the existence of a new type of process that 
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contributes to the quark-antiquark thermal decay width besides the Landau damping. 

We will deal with static quarks only. The static case is relevant also for the study of 
bound states made of quarks with a large but finite mass m, like quarkonia, in a thermal 
medium. Quarkonium is expected to exist in the medium if the temperature and the other 
thermodynamical scales are much lower than m. In this situation, one may consistently 
integrate out the mass from QCD and expand order by order in 1/m. The leading order of 
the expansion corresponds to QCD with a static quark and a static antiquark. Higher-order 
corrections in 1/m may be systematically included in the framework of non-relativistic QCD 
(NRQCD) 

Bound states at finite temperature are systems characterized by many energy scales. 
There are the thermodynamical scales that describe the motion of the particles in the thermal 
bath: the temperature scale T (we will not distinguish between T and multiple of vrT), the 
Debye mass m^), which is the scale of the screening of the chromoelectric interactions, and 
lower energy scales. In the weak-coupling regime, which we will assume throughout this 
work, one has ~ gT ^ T. Moreover, there are the scales typical of the bound state. 
In the case of a system of two static sources, the scales may be identified with the inverse 
of the quark-antiquark distance r and the static potential, the first being much larger than 
the second. We will also assume that 1/r is much larger than the typical hadronic scale 
Aqcd; i-e. we will concentrate on the short- distance part of the potential. This may be 
the part of the potential relevant for the lowest quarkonium resonances like the J/ip or the 
T(IS'), which are the most tightly bound states. Thermodynamical and bound-state scales 
get entangled and different hierarchies are possible. Bound states are expected to dissolve in 
the bath at temperatures such that m^j is larger than the typical inverse size of the bound 
state. Hence we will concentrate on the situation 1/r > mo (i-e. 1/r ^ mo or 1/r ~ mo), 
and distinguish between the two cases 1/r ^ T and 1/r ^ T. 

The paper is organized as follows. Sections [III and IIIII are introductory: they deal with 
QCD with static sources, which we call static QCD for short, but do not include bound 
states. In Sec. HIl we write the quark and gluon propagators in static QCD at finite T. In 
Sec. Illlt we summarize one-loop finite T contributions in static QCD that are relevant to the 
present work. In Sec. lIVl we introduce the relevant EFTs and calculate the static potential in 
a situation where the inverse distance between the static quark and antiquark is larger than 
the temperature of the thermal bath: 1/r ^ T ^ m/j. When T is as small as the binding 
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energy, we also calculate the leading thermal contribution to the static energy and the decay 
width. In Sec. |Vl we provide an alternative derivation of the potential in perturbative QCD. 
In Sec. rvT| we calculate the static potential in the situation T ^ 1/r > m^i. In Sec. IVIIt 
we summarize and discuss our results and list some possible developments. 



II. STATIC QCD AT FINITE T 

We consider here QCD with a static quark and antiquark; in particular, we write the 
quark, antiquark and the gluon propagators at finite T. To simplify the notation, we will 
drop the color indices from the propagators. Throughout the paper, the complex time 
contour for the evaluation of the real-time thermal expectation values goes from a real 
initial time ti to a real final time tf, from tf to tf — iO^, from tf — iO^ to — iO^ and from 
ti — iO^ to ti — i/T. The propagators will be given with this conventional choice of contour. 
Furthermore, the following notations will be used. We indicate thermal averages as 

(o) -Z^i^:!^ m 

where H is the Hamiltonian of the system. We also define 

Mk') = ^g^^, (2) 
Mk') = (3) 



A. Quark propagator 

In order to show the behaviour of static sources in a thermal bath, it may be useful to 
consider first a quark (or antiquark) with a large but finite mass m, m ^ T, and then 
perform the m — oo limit. 

We define the propagators 

s^x) = {Mx)4mT, (4) 

S<^ix) = -{^ljliO)Mx))T, (5) 

where ip is the Pauli spinor field that annihilates the fermion (in the following, the Pauli 
spinor field that creates the antifermion will be denoted x)- The free propagators, 

C>(0) _ ^ ^C>(0) q<{0) _ ^ ^C<(0) (f^\ 
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A;05>(°)(A;) =m^>(°)(A;), (7) 
k'S<^'^\k)=mS<^'^\k), (8) 

where we have neglected corrections of order 1/m or smaller: they will eventually vanish in 
the m ^ oo limit. 

If the heavy quarks are part of the thermal bath, they satisfy the Kubo-Martin-Schwinger 
relation: 

5<(0)(fc) = _e-fcVT^>(0)(^), (9) 

From the equal-time canonical commutation relation it follows the sum rule 

— (S>^'\k)-S<^'\k)) = l. (10) 

The solutions of the equations fl7|)- f|T0|) are 

^>(o)(fc) = (l-nF(A;°))27r5(fc°-m), (11) 
5<(o)(;t) = -nF(A;°)27r(5(A;°-m). (12) 

The spectral density pp •* is given by 

pP{k) = S>^°\k) - S<(°)(^) = 27r(5(fc° -m), (13) 

and the free propagator, 

= - ^(-x°)5<(°)(a;), (14) 

is given in momentum space by 

S^'\k) = ^^^^ - nAk') 2n5ik' - m). (15) 

In the static limit m oo, the propagators simplify because nF{m) — for m — * oo. 
Moreover, we may get rid of the explicit mass dependence by means of the field redefinition 
ip 'i/'e"*™*, which amounts to change k^ — m to k^ in the expressions for the propagators 
and the spectral density; they read now 

5>(°)(A;) = 27r5(fc°), (16) 



le 



pW(A;) = 27r5(A;°). 



(17) 
(18) 
(19) 



The free static propagator is the same as at zero temperature. On the other hand, if we 
would have assumed from the beffinnine; that S<^^\k) = 0, i.e. that there is no backward 
propagation of a static quark (in agreement with the Kubo-Martin-Schwinger formula in the 
m ^ oo limit) then, together with the equations of motion k'^S^^^^k) = 0, k'^S'^^^^k) = 
(obtained after removing m via field redefinitions) and the sum rule (fTOl) . we would have 
obtained Eqs. ([TB]). ([T5D and (1191). 

The real-time free static propagator for the quark reads 



and for the antiquark 



^W(A;) S<W(A;) 
S>io){k) {S^^\k)y 

( ' 

v 



Jal3 



P + ie 
2-n5{k^) 





—i 



\ 



(20) 



k^-iej 



sS(fc) = 5, 



al3 



2Tx5{k^) 





—i 



-k^ 



(21) 



The main observation here is that, since the [S^'j](A;)]i2 component vanishes, the static quark 
(antiquark) fields labeled "2" never enter in any physical amplitude, i.e. any amplitude 
that has the physical fields, labeled "1", as initial and final states. Hence, when considering 
physical amplitudes, the static fields "2" decouple and may be ignored. 
The propagator s[^°^ may be written in a diagonal form as 



k^ + ie 




U(o) 



(22) 



where 



U(o) 




and for further use [U^'^-'l ^ 




(23) 



Throughout the paper, we will use bold-face letters to indicate 2x2 matrices in the real-time 
formalism. 
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B. Gluon propagator 

The gluon propagator in the real-time formalism can be written as 



where 



D>{k)= / d^xe'''-^{A^{x)Ap{0)) 



T, 



Dc.p{k) = [ d'^xe"'-^ [e{xo){A^{x)Ap{0))T + 9{-xo){Ap{0)Ao.{x)h 



Gluons being bosonic fields, the Kubo-Martin-Schwinger relation reads 

D<ik) = e-^°/^D>ik), 



from which it follows that 



where 



D>p{k) = {l + nB{k°))pB^p{k), 
D<f,{k) = nB{k'')pBafs{k), 

PBapik) = D>^ik) - D<^ik), 



is the spectral density. 

We may express D^^ also in terms of the retarded and advanced propagators D 



D%{k) = - I d'xe''-^e{-x,){[A^{x),A0)])r. 



we have 



PBo.m=Dlp{k)-Dt^{k), 

Dap{k) = D%{k) + D<^{k) = Dt^{k) + Dl^{k) 



^ + {^ + nB{k ) ) PBa/3(A;). 



(35) 



In the free case, in Coulomb gauge, the longitudinal and transverse propagators have the 
following expressions 







dS?(A:) 



A;2 

-I. 



(36) 



A;2 



A;2 + ze 
^(A;0) 27r5(A;2) 



^(-A;0) 27r5(A;2^^ 



A;2 



+ 27r5(A;2)nB(|A;°|) 



1 1 
1 1 



(37) 



where k"^ = (k^)'^ — k"^ . Note that the longitudinal part of the gluon propagator in Coulomb 
gauge does not depend on the temperature. The temperature enters only the transverse 
part, which splits in the sum of a T = piece and a thermal one. 



C. Lagrangian 



The Lagrangian of QCD with a static quark, a static antiquark and nj massless quark 
fields qi is 

Uf 

C = -^F^^F^^-^ + 5^ tpq, + ^HDo^ + x^tD^x, (38) 
1=1 

where iD'^ = ido — gA'^, iD = iV + gA, and igFfj^^ = [D^,D,^]. The free static quark 
propagator is given by Eq. fl2U]) . the free static antiquark propagator by Eq. fl^ and the 
free gluon propagator (in Coulomb gauge) by Eqs. (|36|1 and (l37|l . Note that transverse 
gluons do not couple directly to static quarks. 



III. ONE-LOOP FINITE T CONTRIBUTIONS IN STATIC QCD 

Throughout this work, we will assume that T, gT 3> Aqcd! this enables us to evaluate 
thermal properties of QCD in the weak-coupling regime. In this section, we consider one- 
loop thermal contributions to the static quark propagator, quark-gluon vertices and gluon 
propagator. When the loop momenta and energies are taken at the scale T and the external 
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momenta are much lower, so that we may expand with respect to them, these correspond 
to the hard thermal loop (HTL) contributions jl3| . 



FIG. 1: Self-energy diagram. The continuous line stands for a static quark, the dashed one for a 
longitudinal gluon. 



FIG. 2: Vertex correction to the static quark line; the incoming gluon is longitudinal. 




FIG. 3: Vertex correction to the static quark line; the incoming gluon is transverse. 

The one-loop contributions to the static-quark self energy, to the static-quark lon- 
gitudinal-gluon vertex and to the static-quark transverse-gluon vertex are displayed in 
Figs. [H [2] and [3] respectively. It is convenient to fix the Coulomb gauge. In that gauge, lon- 
gitudinal gluons do not depend on the temperature (see Eq. (!36l) ) and the above diagrams 
do not give thermal contributions. Moreover, if evaluated in dimensional regularization they 
vanish after expanding in the external momenta. Throughout this work, we will adopt the 
Coulomb gauge unless stated otherwise. 

Momenta and energies of order T contribute to the gluon self-energy diagrams. Since 
only longitudinal gluons couple to static quarks, we will focus on the longitudinal part of 
the polarization tensor. This will be the only component of the gluon polarization tensor 
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relevant to the paper. Diagrams contributing to the thermal part of the longitudinal gluon 
polarization tensor in real-time formalism at one-loop order are shown in Fig. |H In Sec. 1111 Al 
we will give a general expression, and in Sees. 1111 Bl 1111 C\ 1111 Dl we will expand it in some 
relevant limits. 



FIG. 4: Diagrams contributing to the longitudinal component of the gluon polarization tensor 
at one-loop order. The continuous loop stands for light (massless) quark loops, dashed lines for 
longitudinal gluons and curly lines for transverse gluons. Ghosts do not contribute to the thermal 
part of the gluon polarization tensor 



I 
I 
I 
I 

I 




A. The longitudinal gluon polarization tensor 



Summing up all thermal contributions from the diagrams of Fig. HJ we obtain (for details 
see Q): 



[Km 



thermal 



thermal 



noo(A;° + «e, k) 
noo(fc°-^e,^) 



thermal 



thermal 



(39) 
(40) 



n00(fc)]thermal " ^00, F(fc)]thermal + PoO, G (fc)] thermal ' 
„2t^ „ r+co 



(41) 



00, F 

thermal 



X 



27r2 
4(g")^ + k 



dq'\q'\np{\q'\) 



— CO 

0\2 I l2 a^OuO 



+ 



4|g0||fc| 

4(g°)2 + A;2 + 4gQfc° 
4|g0||fc| 



In 



4gOA;0 \ - 2q^ + 2\qyk\ 

P_2g0fc0_2|g0||^| 

P + 2gOA;0-2|gO||fc| 



In 



A;2 + 2gOA;0 + 2|gO||A;| 



(42) 
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thermal 



1 + 



27r2 

(2g° - 
8(g0)2 



dq'\q>Bi\q'\) 



1 
2 



2(g' 



o^2 



\k\ i\k\^ + {q 



0\2\2 



2|g0| 



8|gO|3|A;| 



X In 



(2g0 


- A;0)2 


4(g0 


- A;0)2 


1^1- 






1^1 + 







0\2n2 



(2g0 - 



1 



2|gO||A;| 



X In 



A;2 - 2q^k^ 



8|gO|3|A;| 



2q^ky ^ e 




2gOA;0 



8|gO|3|A;| 

2|g°||^| 



2|gO||A;| 




(43) 



p_2gOA;0_2|gO||^| 

where "R" stands for retarded, "A" for advanced, "F" labels the contribution coming from 
the loops of Uf massless quarks (first diagram of Fig. Hj) and "G" labels the contribution 
from the second, third and fourth diagram of Fig. HI A'c = 3 is the number of colors and 
Tp = 1/2. In the context of the imaginary-time formalism, Eas. f|^2l) and (H3l) can be found 
also in textbooks like 15|. The original derivation of ( H3l) is in |l6 |. 

The retarded and advanced gluon self energies contribute to the retarded and advanced 
gluon propagators. From the retarded and advanced gluon propagators we may derive the 
full propagator, the spectral density and finally all components of the 2x2 matrix of the real- 
time gluon propagator along the lines of Sec. IIIBI In the following, we study Eqs. ( j39l) - (l43ll 
in different kinematical limits. 



B. The longitudinal gluon polarization tensor for <^T 



The typical loop momentum is of order T. If we expand [n^Q(A;)]^j^^^^^j and 
in k^ <^ T \k\ and keep terms up to order the result is 



thermal 

Re Wo(A;)]^„ = Re [n^o(fc)] thermal = 



rdq'q'nAq') 
^ Jo 



^ Jo 



dq\^n^{q^) 



1 - 



\k\ - 2g° 



\k\ + 2q^ 



\k\ 



In 



\k\ - 2g° 



|A;| +2g0 



(44) 
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Im [Ulik)] 



thermal thermal 

^ k J\k\n 



Ifcl 



?^b(|^|/2) 



dq'q'uBiq' 



\k\/2 



.(45) 



Equation ( l45l) and the gluonic part of ( H4|) are in agreement with 



ith Q. 



C. The longitudinal gluon polarization tensor for ~ ^ T 

If we assume that all components of the external four-momentum are much smaller than 

the loop momentum g° ~ T, then we may expand [noo(^)] thermal [■'^oo(^)] thermal ~ 
\k\ <^ T. At leading order, we obtain the well-known HTL expression for the longitudinal 

CI 

gluon polarization tensor, which may be found, for instance, in [9|: 



k' 



p - \k\ 



(46) 



I- Wo(A:)] thermal = (^)] thermal = ^ ^ ^(-^^^ ' ^^^^ 

where m^) is the Debye mass: 

ml = ^iN, + TFnf) . (48) 

We have used that / dq^ q^ n^iq^) = and / dq^ q^ n^iq^) = n^T^/Q. Note that 

Jp Jo 
the expansions for \k\ — of (H4|) and (H5|) and those for — of fj46l) and (H71) agree with 

each other at leading order. 

i. The resummed longitudinal gluon propagator 

The longitudinal polarization tensor induces corrections to the longitudinal gluon prop- 
agator: 

D^\k) = ^ - ^nSj^(fc) + . . . . (49) 

Since H^q^ contains a real and an imaginary part, also -D^'^ acquires a real and an imaginary 
part. 



12 



If the typical momentum transfer is of the order of the Debye mass, \k\ ~ m^,, then the 
series in (l49ll needs to be resummed: 

D^\k) = ^ (50) 

The resummed longitudinal propagator depends on and has a real and an imaginary 
part. The Debye mass rriu plays the role of a screening mass for longitudinal gluons whose 
momenta are such that <^T and ~ mp . A study of the resummed gluon propagator 



in the real-time formalism may be found in 



17| 



The role of the screening mass can be made more evident if we assume further that 
/c°<C|fc|~m£)<^T and expand Eq. fISUl) in k^ up to order then we obtain 

D^^{k) = ^ ±-t ^^^^ , (51) 

k^ + ml 2 \k\ (^p + ^2^j 

where the "+" and "— " signs refer to the retarded and advanced propagators respectively. 
The corresponding spectral density is 

PBoo(A:) = D^k) - D^,{k) = vr ^ , • (52) 

Following Sec. IIIBI and expanding in k^ also the Bose factor, nB{k^) ~ T/k^, at leading 
order in A;°, we obtain that the resummed HTL longitudinal propagator in the real-time 
formalism is 





k^ + mi yo -ij \k\ (p , _2 

D. The longitudinal gluon polarization tensor for ^ T ^ 

If we assume that \k\ ^ T ^ k^, then the expression for the longitudinal gluon polariza- 
tion tensor may extracted from Eqs. (jH]) and fH5l) by expanding for large \k\/T. At leading 
order, we obtain 

[noo(^)] thermal = [^oo (^)] thermal = ]^ " (^^) 

The result is real and does not depend on k. Moreover, only the gluonic part of the polar- 
ization tensor contributes in this limit and at this order. Higher-order real corrections are 
suppressed by T^//c^, while higher-order imaginary corrections are exponentially suppressed. 
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IV. BOUND STATES FOR 1/r > T 



Starting from this section, we shall address bound states made of a static quark and 
antiquark in QCD at finite T. Bound states introduce extra scales in the dynamics, besides 
T and m^), that we have to account for. The most relevant one is the distance r between 
the quark and the antiquark. Throughout the paper, we will assume that 1/r ^ Aqcd- We 
will further assume that also the binding energy of the quark-antiquark static pair is larger 
than Aqcd- 

First, we deal with the situation where the inverse distance of the two static sources is 
much larger than the temperature: 1/r ^ T. Under this condition, we may integrate out 1/r 
from static QCD at T = order by order in Og. The EFT that we obtain is potential non- 
relativistic QCD (pNRQCD) in the static limit 
as 
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19| . whose Lagrangian can be written 



4 



i=l 



(fr Tr < 



i9n + C^'^ 



S + O^ 



1 «Vo 

2N^ r 



O 



+l^ATr|0V-^iS + SV-^E0| + ^Tr|o¥-^EO + O^Of-c/E| + ... . (55) 

The fields S = S %/ \/Wc and O = 0°^ T"/ \/Tp are static quark-antiquark singlet and octet 
fields respectively, E is the chromoelectric field: i?* = and Cp = {N^ — l)/{2Nc) = 4/3. 
The trace is over the color indeces. The matching coefficients ay^, Va, Vb are at leading 
order: ay^ = as, ay^ = Og, Va = 1 and Vb = 1- Gluon fields are multipole expanded and 
depend only on the center of mass coordinate; they scale with the low-energy scales (T, rri^), 
Os/r, Aqcd; •■■) that are still dynamical in the EFT. The dots in the last line stand for 
higher-order terms in the multipole expansion. 



A. Singlet and octet propagators 



The free real-time singlet and octet static propagators at finite T are similar to the free 
static quark propagator fl20|) . although singlet and octet are bosons: 



gsinglet^^) 



p° + Cpayjr + ie 
I 27r5(/ + Ci.ay7r) ' -) 





—i 



\ 
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U(o) 



p° + CpOivJr + it 




U(o) 



(56) 



p 



'° + CpOLyJr - iej 





—i 



\ 



-avJ{2Ncr)+ie 
/ 



-avJ{2N,r)+te 




— « 



u 



(0) 



(57) 



\ p''-avJ{2N,r)-ie/ 
B. Non-thermal part of the singlet static potential 

The contribution to the singlet static potential coming from the scale 1/r can be read 
from the Lagrangian (!55l) : it is just the Coulomb potential, which in real-time formalism 
reads 

ay^{l/r) I 1 
-1 



Vs(r) = -Cf- 



(58) 



where ay^(l/r) is a series in as', at leading order, ay^(l/r) = as{l/r). Starting from order 
Og, ay is infrared divergent; these divergences, which appear at zero temperature, have 



been considered elsewhere 



20| and will not matter here. 



The matrix in (135]) is such that 







(i o\ 


[o -I) 




^0 -l) 



[U 



(0)1-1 



(59) 



Together with Eq. fl22l) . this guarantees that 

\ 



p^ + ie 
2vr5(pO) - 
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71=0 



-«Vs(r)) 



p^ + ie 
27r5(/) 



\ 

—i 



p^ — It 



(60) 



p^ - It, 

i.e. that, like at T = 0, the sum of all insertions of a potential exchange between a free 
quark and antiquark amounts to the propagator fISB]) . 

The singlet static potential does not get only contributions from the scale 1/r, but, if the 
next relevant scales are the thermal ones, it will also get thermal contributions. These will 
be the subject of the following sections. 
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C. Chromoelectric correlator 



In the paper, it will become necessary to calculate the chromoelectric correlator 
{E°'{t)(j){t,0)'^'l^E''{0))T where 0(^,0)^^-' is a Wilson line in the adjoint representation con- 
necting the points {t, 0) and (0, 0) by a straight line. Such a correlator enters each time we 
consider diagrams with two chromoelectric dipole insertions. In the real-time formalism, the 
chromoelectric correlator is a 2 x 2 matrix in the field indices "1" and "2" . We shall write 
it as 

where, at zeroth order in as, Df^u{k) is the free gluon propagator. In Coulomb gauge, the free 
gluon propagator has been given in Eqs. ( !36l) and ( !37l) : since the chromoelectric correlator 
is a gauge invariant quantity the choice of the gauge does not matter. At zeroth order in 
as, the thermal part of {E''{t)(f){t,0)'fjE\0))T is 



(62) 

Note that for t = Eq. ( l62l) gives the thermal part of the gluon condensate in the weak- 
coupling regime: 

OO 

3^ _ _4t-t4 



d^k . / 1 1 



where we have used / dk k ub 







Boltzmann law (see, for instance, [2l|). 



k) = TT T /1 5. Equation (1^ agrees with the Stefan- 



D. Thermal corrections to the singlet static potential 

We calculate now the leading thermal contributions to the static potential assuming for 
definiteness that T and rriD are the next relevant scales after 1/r, i.e. that the binding 
energy is much smaller than mo- We will further assume that all other thermodynamical 
scales are much smaller than the binding energy, so that we can ignore them. We recall that 
from an EFT point of view, only energy scales larger than the binding energy contribute 
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to the potential, which is the matching coefficient entering the Schrodinger equation of the 
bound state, while all energy scales contribute to physical observables like, for instance, the 
static energy [l^. We will comment on the impact of degrees of freedom with energies and 
momenta of the order of the binding energy in Sees. IIVEI and lIVFi 

The calculation proceeds in two steps. The ffist step will be performed in Sec. IIVD 1[ 
It consists in integrating out from the pNRQCD Lagrangian (155|) modes of energy and 
momentum of order T. This modifies pNRQCD into a new EFT where only modes with 
energy and momentum lower than the temperature are dynamical. For our purposes, it only 
matters that the pNRQCD Lagrangian gets additional contributions in the singlet and in 
the Yang-Mills sectors. In the first sector, it is an additional contribution to the singlet 
static potential. In the second one, the additional contribution corresponds to the HTL 
Lagrangian ^2^ . 




The second step, which will be performed in Sec. IIVD 2t consists in integrating out from 
the previous EFT modes of energy and momentum of order niD. The resulting EFT will 
only have degrees of freedom that are dynamical at energy and momentum scales lower 
than the Debye mass. In the singlet sector, the EFT gets modified by a further additional 
contribution to the static potential. 

In summary, if both T and m^) are much larger than the binding energy, the real-time 
Coulomb potential fl58|) . gets two type of corrections 5Vs(r). The first one comes from the 
scale T and the other one from the scale rriD- 

1. Contributions from the scale T 

The leading thermal correction 5Vs(r) is induced by the dipole terms O^r ■ gE^ + S^r*- 
gEO in the Lagrangian ( |55l) . It reads (see 13] for the T = case): 





[5V,(r)],2 = 0, 

2 

[SY,{r)],, = ^g'CF^^l■ 



dt e 



-itAV 



mmt,o):fE\o))T 



(65) 
(66) 



J 11 



d'^k 



27i5{-k°-AV) ik^f'Du{k) + pBoo{k) ,(67) 
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where AV = - ( — - + Cpay ] ~ — ^— ^- The corresponding Feynman diagram is shown 
r \2Nc / 2r 

in Fig. [5l Integrals over momenta have been regularized in dimensional regularization {d is 
the number of dimensions, fi is the compensating scale). In Eq. flMl) . i/{—k^ — AV + ie) 
is the "11" component of the static octet propagator; Eq. (!66|) vanishes because the "12" 
component of the static octet propagator vanishes and in Eq. fl^7|) . 27[6{—k^ — AV) is the 
"21" component of the static octet propagator. Note that vertices of type "1" and "2" 
have opposite signs. Equation (!65l) . which may also be read [—i6Y s{r)]22 = [~'^<^'^s('")]ii, 
reflects the relation existing between the "11" and "22" components of the propagators in 
the real-time formalism. 




FIG. 5: The single continuous line stands for a singlet propagator, the double line for an octet 
propagator, the circle with a cross for a chromoelectric dipole vertex and the curly line connecting 
the two circles with a cross for a chromoelectric correlator. 



We are interested in calculating the contribution to the integrals in Eqs. ( !64l) -( !67|) from 
momenta /c ~ T. Since T ^ AV, we may expand in AV/T. Moreover, at leading order, 
the propagators in Eqs. flMl) and (jnZD are the free ones, Dqq^ and D.°\ given in Eqs. (15^ 
and fl37j) . However the leading-order thermal contribution, which would be of order 
vanishes: 

[S^s{r)l, = -z/ ^f^'-' J _^oV ^e ^^^^^'^ ^^^1^°') = ° ' ^^^^ 
[6Y^{r)],, = ./C^^/i^-^ (05 2^'^(-^°) 47r5(fc^)nB(|A:°|) = 0. (69) 

Several next-to-leading order corrections are possible, because several scales are still dy- 
namical in the EFT: we may have corrections of relative order AV/T, rrin/T, (rT), ag and 
so on. 

(1) First, we consider corrections of order AV/k^ or higher to the quark-antiquark prop- 
agator, which contribute to order g'^r^T^ x AV/T or higher to 5Vs(r): 

l*V,(.)|„ ^ -./C.^,-/ ^^-^,-L_ [(,0)=D«»'(.) +PdS(*)]^^ 
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= ^Cp^r'T' AV fiAV/T) - i\ a, {AVf n^{AV) , (70) 

= Cf tts {AV fn^{AV) = -2 ^ Im [5V,(r)],, , (71) 

where 

f{z)= / rfx^P^ (72) 

and P stands for the principal value. Since T ^ AV, we can expand the above expressions 
in AV/T obtaining 

5Y,{r) = ^-N,CFalrT^ i^^ + . . . - '-NlCp alT ^ j + " " " " (^3) 

The dots stand for higher-order real and imaginary terms. We note that the matrices in 
Eq. (!73l) are such to be diagonalized by the matrix [U'-'^^]^^: the first one is the same as in 
f l59ll . the second one satisfies 

= |U<°>|-' I ^ "1 |U<°>|-'. (74) 

V" V 

The leading real contribution in Eq. (1731) is of order g'^r'^T^ x AV/T and comes from 
the diagram shown in Fig. [61 The leading imaginary contribution in Eq. (!73l) is of order 
g2^2rp3 ^ [AV/T)"^ and comes from the diagram shown in Fig. [71 This imaginary part 
is a very peculiar feature of QCD. It originates from the fact that thermal fluctuations of 
the medium at short distances may destroy a color-singlet bound state in an octet quark- 
antiquark state and gluons. This process is kinematically not allowed at zero temperature 
where only static octets may decay into singlets with a leading order width F = a^/ (12r). 
In some kinematical situations, the role of singlet-octet transitions in quarkonium-hadron 




scattering in strongly interacting matter has been discussed in Ref. [23 1. 

Here we have assumed T ^ AV. In the kinematical situation T ~ AV, Eq. f[70|) would 
provide a contribution to the static energy rather than to the static potential of the quark- 
antiquark pair. This contribution will be discussed in Sec. IIVEI 

(2) Another source of next-to-leading order corrections comes from next-to-leading order 
corrections to the chromoelectric field correlator. Their contribution to the static potential 



at zero temperature has been considered in 2J] . They also contribute at order g"^ x a, 
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FIG. 6: Feynman diagram giving the leading order real contribution to Eq. (|7U|) . The symbols are 
like in Fig. [5l The cross stands for a AV insertion in the octet propagator. 




FIG. 7: Feynman diagram giving the leading order imaginary contribution to Eq. (jTOp . The 
symbols are like in Fig. [6j 



g2 r'^T^ X [m£,/T)'^ to the thermal part of the static potential. The chromoelectric correlator 
enters in the potential in the expression 



{2'kY -k^ + it 



(A;°)2d,,(A;) + A;^Doo(A;) 



11 



Since i/{-k^ + = -iF{l/k^) +7c6{-k^) and (A;°)2 Dii(fc) + P Doo(A;) is even in k°, 
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only the 7c6{—k^) component of the static quark-antiquark propagator contributes, therefore 



only the limit for k^ ^ oi 



matters. In order to evaluate it, it is 



{ky-Du{k) + k^-Doo{k) 
convenient to perform the calculation first in temporal-axial gauge ^4° = 0. In this gauge, 
the chromoelectric field is simply E = —OqA. Hence all corrections to the chromoelectric 
correlator are encoded in the spatial part of the gluon propagator alone: at one loop the 
correction is provided entirely by the gluon self energy. In temporal-axial gauge, from the 
transversality relation of the polarization tensor it follows that (compare with the expressions 
of the propagators in [3]): 

P 



lim {kyDf,'''{k)\ _ ^ , 

I temporal— axial gauge rh iiqq \^rv j 



lim i- 



(75) 



temporal— axial gauge 



Since in the k^ ^ limit n^Q^(fc) is equal in Coulomb and temporal-axial gauge, we can 
also write that (compare with Eq. (1501) ) 



lim (k' 



0^2^R,A 



(fc) 



lim pD^A^(k) 

fcO^O 



(76) 



Coulomb gauge 



temporal— axial gauge 

The left-hand side is the only term of the chromoelectric correlator contributing to the 
potential in temporal-axial gauge: it may be evaluated by calculating the right-hand side in 
Coulomb gauge. At one loop, the right-hand side gets contribution from the gluon self-energy 
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diagram shown in Fig. [HI hence, at next-to-leading order we can write 



r 



A-d 



d-l 

2 



r 



d-l 



fx 



d'^k 



TTC 



Si-k')k' [SBooik)],, 



27r5(-A;°) k^ [SBooik)] 
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where 



[6Boo{k)],, = (1 + n^{k')) {6D^,{k) - 6D^,{k)) , 



6D^^{k) 



A;4 



(77) 
(78) 

(79) 
(80) 
(81) 



and, the relevant limit for the gluon polarization Il^Q^{k) in Coulomb gauge is given by 
Eqs. dH]) and fH5|) . Finally, the correction to the real-time potential reads 



-|C(3) Cp^r'Tml + |c(3) N^Cp al 








47r 



T2 2 

+ 7i5 + ln7r-ln — H 4 In 2-2 



cm 



9 




(82) 



where e = (4 — c/)/2, 'Je is the Euler gamma and ( the Riemann zeta function (C(2) = v^/6). 
Note that in Eq. fl82p . besides terms that are proportional to the Debye mass there are finite 
terms, both in the real and in the imaginary part, that do not depend on it. 



FIG. 8: The symbols are like in Fig. [5l The dashed blob stands for a one- loop self-energy insertion 
in the gluon propagator. 

Equation flS^ contains an imaginary contribution. The origin of this contribution is 
different from the one in Eq. (175]) . The one here comes from the imaginary part in the 
gluon self energy, which is due to to the scattering of particles with momenta of order T in 
the thermal bath with space-like gluons, (/c°)^ < (Landau damping) while the one in 
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Eq. (173!) signals the thermal break up of a quark- ant iquark color singlet pair into an octet 
one. 

The result in Eq. fl82p is infrared divergent and, in an EFT language, calls for an opposite 
ultraviolet divergence from lower energy contributions: the two divergences should cancel in 
all physical observables. In the following section, we will show that the infrared divergence 
generated by the diagram in Fig. [8] when integrated over momenta of order T is canceled 
by an ultraviolet divergence in the same diagram when integrated over momenta of order 
m£,. For the purpose of the cancellation and the calculation of the static potential in 
the situation 1/r ^ T ^ rriD ^ AV^ it is irrelevant how we regularize both divergent 
contributions, as long as they are regularized in the same way: the divergences as well as 
any scheme dependence cancel in the sum. However, it may be that expression ([H2D is used 
in intermediate calculations that require further regularizations. An example is the solution 
of the Schrodinger equation for quarkonium in a thermal medium when the kinetic energy 
of the bound state is larger than or of the same order as rriD so that the cancellation of the 
divergences does not occur at the level of the potential. In such a situation, it is important 
to provide divergent contributions in a standard regularization scheme. Equation ( l82l) has 
been obtained by regularizing in dimensional regularization only the integral in k. This is 
not sufficient to qualify Eq. (182|) as a standard dimensional regularization scheme result, 
which would also require the calculation of the gluon polarization tensor at order e in the 
dimensional regularization expansion. This order, combined with the 1/e divergence of the 
integral in k, contributes to the finite part of ( l82l) . To our knowledge the expression of 
the gluon polarization tensor at order e is not known in the literature and its calculation is 
beyond the purposes of this work. However, it may become necessary for a proper calculation 
of the quarkonium potential at short distances when the Debye mass is of the same order 
as or smaller than the kinetic energy. 

(3) Finally, a source of higher-order contributions comes from higher-order terms in the 
multipole expansion. These are of two types. First, they may involve operators of higher 
order in the multipole expansion in the pNRQCD Lagrangian, like, for example, the operator 
— Tr jOVV-'r'^ gD^D^E^ S}. However, the contribution of this operator vanishes for reasons 
similar to the ones that led to the vanishing of the leading-order thermal contribution in 
Eqs. (1681) and fl69l) . Second, they may involve diagrams with three or more insertions of the 
operators Tr|o^r-(7^s| or Tr |o^r ■ (7^o|. They contribute to order g^r'^T^ x g"^ rT 
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or higher to the thermal part of the static potential and hence are suppressed by at least 
a factor rT with respect to the corrections considered in the previous paragraph. We will 
neglect them in the following. 



2. Contributions from the scale mo 



In the previous section, having integrated out T has led to a new EFT, which in the 
Yang-Mills sector coincides with the HTL EFT and in the singlet sector shows a potential 
that is the sum of the terms (jSHD, (1751) and fl5^ . In the weak-coupling regime, the Debye 
mass m£), which is given by Eq. (HSj) . is smaller than the temperature. We assume that m£, 
is the most relevant scale in the new EFT. The leading contributions to the static potential 
from the scale mr, originate from the diagram shown in Fig.[5]when integrated over momenta 
k mo- They are given by 

d'^k i 



r 



A-d 



d-l 

2 



/i 



A-d 



d-l 
[5V,(r)],, = 0, 



/i 



{2TxY-k^ + ie 
d'^k 



(27r)' 



{kyT>,,{k) + eT>o^{k) 
(A;°)2d,,(A;) + PDoo(A;) 
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d 



4~d 



d'^k 
{2nY 



27rSi~k') iky-Duik) + eBooik) 
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(83) 
(84) 
(85) 
(86) 



where T)fj_^{k) is now the HTL resummed propagator. After integration in k^ only Doo(0, k) 
contributes, the expression of which can be found in (135]) . Substituting and performing the 
dimensional integrals, we obtain 



<5V,(r) 



F 2 3 

— tts r mjj 
6 \ 



1 



i T rrij-, 7^; + In vr + In 

6 




(87) 



Equation (187|) shows that the scale mr, contributes at order g'^r'^mj^ to the real part of the 
potential and at order g"^ T mjj to the imaginary one. 

The Debye mass effectively plays the role of agluon mass; in this sense, the real part of 



87|) agrees with a result that can be found in 



19|. The imaginary part originates, as the 
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one in Eq. (l82i) . from the imaginary part of the gluon self energy. It shows an ultraviolet 
divergence. This cancels against the infrared divergence of Eq. fl82l) . which, we recall, comes 
from the diagram in Fig. [S] when integrated over momenta of order T. We have already 
commented on this cancellation at the end of the previous section, we will add further 
specifications in Sec. IIVFI 



E. Singlet static energy for T < AV 

If T ~ AV there are no temperature dependent corrections to the potential, but Eq. (170|) 
gives the leading thermal correction to the pole of the color singlet propagator. The real 
part of the pole provides the static energy of a color singlet quark-antiquark pair, whose 
leading thermal part is 

2 (T^ 

6E, = - N,Cf r f{N,a,/{2rT)) , (88) 
o vr 



where the function / has been defined in f[72|) . Minus twice the imaginary part of the pole 
provides the color singlet thermal decay width, whose leading contribution, due to the decay 
of a static quark-antiquark color singlet into a quark-antiquark color octet, is 

n^{N,aJ{2r)) . (89) 



6 r 

In the situation T ^ AV, the thermal width is exponentially suppressed, while the thermal 
contribution to the static energy becomes 

^Es = r3 = -1 vr ^ (E'^(O) ■ E'^(O)). . (90) 



45 3 N, 



thermal part 



Equation (pOj) agrees with the analogous expression for the leading gluon condensate cor 



rection to the quark-antiquark static energy at zero temperature that was derived in 
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F. Summary and comments 

As in the zero temperature case also in a thermal bath, the computation of the real-time 
potential between a static quark-antiquark pair requires integrating out all modes of energy 
and momentum larger than AV . Modes of energy or momentum of order AV or smaller 
enter in physical observables, but, since the y d epend on the binding energy, do not belong 
to a proper EFT definition of the potential [10|. 
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If 1/r ^ T ^ rriD AV and in the weak-coupling regime, the static potential of a 
quark-antiquark pair is obtained by adding Eqs. fl58l) . fl73|) . f l82l) and fl87|) . In its diagonal 
form, it is given by 



The leading term is provided by the Coulomb part: —CfO-s/t. In the real part, the first 
thermal correction is of order g'^r'^T^ x AV/T, the second and third ones are of order 
g2^2rp3 ^ (^rriu/T)'^ and the last one is of order g'^r'^im?j^. The correction proportional to 
asv"^ m?j^ is suppressed by raDjT with respect to the one proportional to agr'^Tm'jj and 
may be neglected. In the imaginary part, the first term is of order g'^r'^T^ x [AV/TY and 
the other ones are of order g"^ r'^T^ x [rriD/T)'^. The dots in (!92l) stand for higher-order real 
and imaginary terms. Temperature dependent higher-order terms are suppressed by powers 
of AV/T, mr)/T, rT, AV/mD and rnio- 

The imaginary part of Vs{r) has two origins. The first term comes from the thermal 
break up of a quark-antiquark color singlet state into a color octet state. The other terms 
come from imaginary contributions to the gluon self energy that may be traced back to 
the Landau damping phenomenon. Both are thermal effects: the first one is specific of the 
non-Abelian nature of QCD, the second one would also show up in QED, although in QED 
the photon polarization tensor would get only fermionic contributions. Having assumed 
T ^ ^ AV^, the term due to the singlet to octet break up is suppressed by {AV/m£,Y 
with respect to the imaginary gluon self-energy contributions, which provide therefore the 
parametrically leading contribution to the imaginary part of the potential. 

The thermal part of Eq. (!92|) is finite because, under the condition l/r^T^m^)^ 
Al^, it provides the leading thermal contribution to the real-time energy and to the decay 




(91) 



where 




^ F 2 S 

+ — as r irq^ + . . . 



(92) 
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width of a static quark and antiquark pair in a color-singlet configuration. Note, however, 
that the non-thermal part of the potential, ay,-, develops infrared divergences starting from 
order ctg, which eventually cancel in physical quantities against contributions from the scale 
Ay. These have been most recently considered in [24]. If 1/r ^ T ^ > m^), the static 
potential is the sum of Eqs. (jSHD, dZSD and (15^ . The thermal part is infrared divergent. Di- 
vergences cancel in physical observables against thermal contributions coming from the scale 
tud- Finally, if 1/r ^ AV^ > T, the static potential is just its Coulomb part fl58l) . Thermal 
contributions affect physical observables through loop corrections that involve momenta or 
energies of the order of the binding energy or smaller, but do not modify the potential. In 
this case, the leading thermal effect on the static energy of the color singlet quark- antiquark 
pair can be read in flSSl) and the leading thermal width, which is due to the singlet to octet 
break-up phenomenon, in fl89|) . 

Our result is also relevant for the case of a quark and antiquark with a large but finite 
mass m. This would correspond to the actual case of heavy quarkonium in a thermal 
medium. However, in the case of finite mass, the relevant scales of the bound state are 
dynamical and the above discussion gets modified accordingly. For a comprehensive review 
in the T = case we refer to [lO]. Bound states of heavy quarks are characterized by the 
energy scales m, mv, mv"^, where v is the relative velocity of the heavy quarks: mv is of the 
order of the inverse size of the bound state and mv'^ is of the order of the kinetic or binding 
energy. The equation of motion of the quark-antiquark bound state reduces in the non- 
relativistic limit to the Schrodinger equation. The potential is the interaction term entering 
the Schrodinger equation. Only modes of energy and momentum larger than mf ^ contribute 
to the potential. With these specifications, Eq. (!92l) provides the static potential of a heavy 
quarkonium in a thermal medium under the condition m ^ mv ^ T ^ rriD ^ mv"^. If 
m ^ mv ^ T ^ mv"^ > mo, the heavy quarkonium static potential is given by the sum 
of Eqs. fl58|) fl73|) and fl82|) only. The potential turns out to be infrared divergent also in 
its thermal part and the comment at the end of paragraph (2) of Sec. IIVD II applies. The 
divergence will eventually cancel in physical quantities against contributions coming from the 
scale mv"^ or lower. If m ^ mv ^ mv"^ > T ^ m^,, the heavy quarkonium static potential 
is just the Coulomb potential (!58l) and thermal corrections enter physical quantities in loop 
involving momenta at the scale mv"^ or lower. 

In the next section, we provide a derivation of Eq. 0921) that does not make use of the 
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EFT language and follows directly from a calculation of the potential in perturbative QCD. 



V. SHORT-DISTANCE THERMAL CORRECTIONS TO THE POTENTIAL IN 
PERTURBATIVE QCD 

In this section, we ask the question of what would be the origin of the thermal part of 
the potential given by Eq. (1921) if we would not introduce any EFT treatment, but simply 
perform a calculation in perturbative QCD under the condition that 1/r ^ T and that the 
binding energy is much smaller than mo- The answer is that the thermal part of fl92|) would 
originate from the longitudinal gluon exchange, with a self-energy insertion, between a static 
quark and a static antiquark shown in Fig. [9] and from the box diagram shown in Fig. [TTl 

1 

I 
I 
I 




I 
I 
I 
I 



FIG. 9: Longitudinal gluon exchange between a static quark and a static antiquark; the dashed 
blob stands for the gluon self energy. 

We first consider the diagram in Fig.[9l which contributes to the physical "11" component 
of the static potential by 

[SV.ir)],, = I -^-^^ e-^'^-^^g' Cp [^5Doo(0, k)] , (93) 

where 5Doo(fc) is defined in Eqs. (1751)- (PT]) and depends on the gluon polarization II^q^. Note 
that we have set to zero the fourth-component of the momentum in the longitudinal gluon: 
corrections would be suppressed by powers of kf/\k\ ~ r or Vs/T or Vs/rrio- Equation 
(!93|) gets contributions from different momentum regions. 

(1) The first momentum region is ~ 1/r. The thermal contribution to the longitudinal 
gluon polarization tensor when |/c|~l/r^Tis provided by Eq. (IMIl . which, substituted 
in Eq. (1^5]) . gives (the integral is finite, hence d = A) 



-^2 



[^Vs(r)], = / -^,e-'\-CF^-^ j = lN^C,alrT\ (94) 
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where we have used that the Fourier transform of 47r/fc'^ is ~r /2. Equation fl9^ agrees with 
the real part of Eq. ( j73l) . 

(2) A second momentum region is \k\ ~ T. Since T -C 1/r, under the condition \k\ ~ T 
we may expand the exponential e'^^'"^ in (l93ll : 



[5V,(r)]i, = /x 



4-d 



(2vr) 



1 



{k ■ f) 



+ 



z5Doo(0,A;) 



11 



(95) 



The first term in the expansion corresponds to a mass correction and cancels against twice 
the thermal contribution of the static quark self energy with a gluon self-energy insertion, 
see Fig. [IDl The second term coincides with the expression in Eq. ( 177|) and gives the same 
result as 



FIG. 10: Gluon self-energy correction to the one-loop self-energy diagram of a static quark. 

(3) Finally, a third momentum region is \k\ ~ mjy. The contribution to the potential 
is like Eq. fl95|) . but now \k\ ~ m£, <^ T and the correct expression for Doo(0, fc) is the 
HTL resummed propagator given in (!53l) . The first term in the expansion corresponds to 
a mass correction, which this time comes from the scale m^i and cancels against twice the 
contribution of the static quark self energy, see Fig. [H when the loop momentum is of order 
rriEi and a HTL resummed gluon propagator is used. The second term gives the same result 

as dHZD. 

We consider now the diagram of Fig. [11], which contributes to the physical "11" component 
of the static potential by (we write the thermal part only) 



(27r)4 ^ ' ' -pO ^i^pO _ qO 

(96) 



The imaginary part of the integral comes from the real part 

Re = Re— ^ = 7r(5(-g°) 
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FIG. 11: Box diagram: the upper line represents a static quark and the lower one a static antiquark. 

which inserted in ( 196|) gives the imaginary part of ( 1731) . The result is exact and does not 
rely on any expansion in the kinematical variables. 

The sum of the contributions coming from the three momentum regions in the integral 
fl93|) and from the imaginary part of fl96|) gives Eq. fl92l) . We note that also the calculation 
in perturbative QCD shows that the imaginary part of Eq. (1731) has a different origin from 
the other ones: it comes from the box diagram of Fig. [Ill which describes a singlet to octet 
to singlet transition, while the other ones come from the gluon self-energy diagram of Fig. [9l 



VI. BOUND STATES FOR 1/r <C T 

In this section, we consider bound states made of a static quark and antiquark in a thermal 
bath at distances such that 1/r ^ T. We still keep that T, 1/r and m^) are perturbative 
scales. We further neglect other thermodynamical scales. 

Under the above condition, the first scale to integrate out from QCD is the temperature 
T. At one loop, this was done in Sec. IIIII Integrating out T leads in the Yang-Mills sector 



to the HTL Lagrangian [2^. In the heavy-quark sector, one loop contributions vanish in the 
static limit. At two loop, see for instance the diagram in Fig. [10], there may be effects. These 
are of order a'^T ~ asT {mo/TY ~ otsma {nici/T) and will be neglected in the following 
where we shall concentrate on the leading contribution coming from the scale mr). 

The next scale to integrate out is 1/r. We assume 1/r > r/iD, and integrate out both 
scales 1/r and m/) at the same time. We shall specialize to the case 1/r ^ m^) in Sec. IVI E[ 
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A. Singlet and octet propagators 



After integrating out the scales 1/r > m^i it may be convenient to introduce quark- 
antiquark fields in analogy with what was done in previous sections. The real-time quark- 
antiquark propagator, S(p), is a 2 x 2 matrix obtained by matching equal time quark and 
antiquark propagators such that [S(p)]jj provides the propagator of a quark-antiquark pair 
of type "z" into a quark-antiquark pair of type "j" . The explicit expressions of the free color 
singlet and color octet quark-antiquark propagators are like those derived in Sec. IIV A[ 



gsinglet (0) ^^-j 



+ ie 
27r5(/) 







-(Octet (0) 



ab 



6, 



ab 



+ it 

27l5{p'') 



p^ — ie) 




^(0) I + 








p^ — ie) 

( ' 

p^ + ie 



p'-' — ie) \ 
Singlet and octet fields have been normalized as in Sec. llVi 

Thermal contributions from the scales 1/r > vhd modify the quark-antiquark propagator. 
In particular, the singlet propagator gets the form 

/ n \ 

gsinglet j^p-j ^ - 5m - V;(r) + 



\ 



V 







u 



(0) 



(97) 



(98) 



— ie) 



\p^ — 6m — Vs{r) + ie p° — 6m* — V*{r) — ie p^ — 6m* — V*{r) — ie) 

^ gsinglet (0) ^ gsinglet (0) • _ ■ y^] gsinglet (0) (p) + . . . ^ (99) 

where in the last line we have expanded with respect to 6m and Vg and introduced the 2x2 
matrices: 



6m 



6m 
—2ilm6m —6m* 



-6m* 



v.. 







-2iimv, -v: I \ -v: 



(100) 
(101) 



B. Matching the mass term 6m 

The static quark (antiquark) self energy at one loop is shown in Fig. [H In the case 
considered here, the loop momentum is of order mn and the HTL resummed gluon propa- 
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gator is used. We match in the real-time formahsm the self energy diagram (normalized in 
color space) with the second term in the expansion fl99l) . S'^™^^'^* (p) [— -j^m] S^™^'''*''°'*(p), 
obtaining: 



A-d 



[5m] 11 = i (ig) CpfJ^' 

= -[(5m]*i , 

[(5m] 12 = 0' 
[5m]2i =2/ C^/i^-'^ 



d'^k 



-k^ + ie -k^-ie 
2n6{-k') [Doo(A;)]ii = 



[Doo(fc)] 



11 



d'^k 



-CpasimD + iT) , (102) 

(103) 
(104) 

[27r5(-A;°) + 27r(5(-fc°)] [Doo(fc)]2i = 2? «s T , (105) 



{2nY 

where i/{—k^ + ie) and ~i/{—k^ — ie) are the "11" components of the static quark and 
antiquark propagators respectively, 2Ti5{—k^) is the "21" component of both static quark 
and antiquark propagators and we have added the contributions from the quark and the 
antiquark; equation (11041) vanishes because the component "12" of the heavy quark and 
antiquark propagators vanishes, see Eqs. fl20|) and fl2T]) . Doo(fc) is the longitudinal HTL 
resummed gluon propagator, whose expression for = is given in Eq. (!53l) . 

The matrix (5m has indeed the form f llOOp : it is diagonalized by the matrix [U'-^-*]^^ and 
5m is given by 

(5m = —Cf tts {jnD + iT) . (106) 

The real part of 5m corresponds to the free energy of two isolated static quarks in the 
imaginary-time formalism, which was first calculated in Ref. 26|] (for further discussions see 
Ref. |27l] ). T he imaginary part of 5m is minus twice the damping rate of an infinitely heavy 



fermion 



C. Matching the singlet static potential Vs 

The matrix elements [Vs]jj are obtained by matching in real time one-gluon exchange 
diagrams that transform a quark- antiquark pair of type "i" into a quark- antiquark pair of 
type "j" with the third term in the expansion ([99]), S'^'^s^''* (p) [-i V^] S^^'^^s^"* (p) . More 
precisely, by matching the diagram of Fig. [12 with [S'''"s'°*(°)(p)] [-i^Au [S''''*^'°*(°)(p)] 
the diagram of Fig. [HI with [S'^'°si^*(°)(p)]22 [-^ Vs]22 [S^'^^sictWj^p^j^^^ ^^le diagrams of 

Fig. [H with [S''*'^si'=*(°)(p)]^J-iVs]i^. [S^''^gi<=t(o)(p)] and the diagrams of Fig. [I5| with 
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^Ssmgiet(o)(^p)j^^ [-iY^].. [S'^i°gi^*(°)(p)]^.^ we obtain 



[Vs(r)] 



11 
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[Vs(r)]i2 = 0, 



[Vs(r)] 
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e-^^-%-g') Cf ( [Doo(0, + [Doo(0, k)]^. 



(27r)3 
-2^Im [V,(r)]n 



(107) 

(108) 
(109) 

(110) 



where the longitudinal HTL resummed gluon propagator, Doo(0, /c), given in Eq. (1531) . comes 
from expanding in the external energy, which is much smaller than the typical momentum 
~ 1/r. 



FIG. 12: Matching condition for [— iVg]]^^; the numbers label the type "1" and "2" propagators. 
All entries in a vertex are of the same type. Vertices of type "1" and "2" have opposite signs. 



[gsinglet (0) (^)] Vs]22 [S^^^S^^* (p)] 



22 



FIG. 13: Matching condition for [— i Vs]22- 



The matching conditions of Figs. fT2] and fT3] fix [Vs(r)]ii and [Vs(r)]22 as in Eqs. (1107^ 
and (11081) respectively. In Figs. [H] and [15] the first two diagrams cancel in the matching 
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FIG. 14: Matching 
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FIG. 15: Matching condition for [— i Vs]2]^. 



against the terms i = j = 1 and z = j = 2 in the sum on the right-hand side. In the case 
of Fig. [lH the last two diagrams vanish, because the "12" component of the static quark 
propagator vanishes, leading to Eq. (11091) : while in the case of Fig. [151 the last two diagrams 
give [Vs(r)]2i as in Eq. ([110]). 

The matrix Vg has indeed the form (11011) : it is diagonalized by the matrix [U*-'^-']^^ and 
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Vs is given by 



ysir) = / 7 — ttC -Ci^^^ l-tCp — m 



r rmD + 1) 

The expression of Vs{r), which we have derived here in real-time formahsm, agrees with the 
analogous expression derived in imaginary-time formalism, after analytical continuation, in 
[sl- It should be emphasized that under the condition 1/r ~ ttid the real part of (11 111) is 
of order agm^), hence sub leading with respect to the imaginary part, which is of order agT: 
the quark-antiquark pair decays before forming the bound state, whose typical time scale is 
proportional to the inverse of the real part of the potential. 

Note that the short-distance expansion of Eq. (II lip would give, up to order r°, the 
Coulomb potential and an r-independent term, Cpo^si'tTT'D + iT), which would cancel the 
mass term derived in (I106p . In the kinematical situation discussed here, this reflects the 
analogous cancellation between the mass correction and the potential correction induced by 
the scale mo discussed in paragraph (3) of Sec. |Vl 

In this section, we have assumed that the singlet propagator has the form (1991) and verified 
that this is indeed the case by performing the matching. We could also have proceeded with 
the reverse logic used in the rest of the paper: by matching 5m and 5Vs we would have 
realized that these matrices are diagonalized by [U^^-']^^, hence, that they fulfill Eq. (!60|) 
and give rise to the resummed propagator fl99l) . 



D. Singlet static energy for 1/r ~ m£) 



Adding the real parts of Eqs. (I106P and (II lip gives the leading static quark-antiquark 
energy for 1/r ~ mo'- 

E, = -Cf a, mD-Cp — e"^"" , (112) 

r 

and the imaginary parts of Eqs. (I106p and (II lip provide the leading static quark-antiquark 
thermal decay width: 



r = 2CFasT 



The thermal width originates from the imaginary part of the gluon self energy. Singlet 
to octet transitions contribute to the decay width as well and the leading contribution is 
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provided by the diagram of Fig. [TTl which gives (see the imaginary part of (173|) ): ST = 
-N^Cf a^T. This contribution is parametrically suppressed by a factor a^, whose natural 
scale is of order 1/r, with respect to the one in (11131) . Note, however, that 6r may be 
numerically as large as 50% of F for as ^ 0.3 and even larger than F for as ~ 0.5; see also 
Fig. [IS 

The static energy given by Eq. (11121) coincides with the leading-order result 27|] of the 
so-called singlet free energy first introduced by Nadkarni J29j (the heavy quark-antiquark 
free energy was defined by McLerran and Svetitsky in |30|) and also studied in lattice 
QCD (see e.g. [3l|, |32|] and [33| for reviews). We recall that the free energy describes a 
thermodynamical property of the system and it is computed from the static quark-antiquark 
propagator evaluated at the imaginary time 1/T (for large temperatures this corresponds to 
small imaginary times), while the static energy studied in this work describes the real-time 
evolution of a quark-antiquark pair and it is computed by evaluating the quark-antiquark 
propagator at infinite real times. The thermal decay width ( 1113^ coincides with the result 
of Ref. 



E. The 1/r » mz5 case 

In the 1/r ^ rriD case, but with still larger than the binding energy, the scales 1/r 
and mD are integrated out in two subsequent matchings. First, the matching at the scale 
1/r can be done in close analogy with the discussion in sections I VI AtTVTCl The potential is 
given by Eqs. (I107l) - (lll0p : since \k\ ~ 1/r ^ rriD we expand Doo(0, k) in powers of ml/k'^. 
We need to regularize the integrals because after expansion they become infrared divergent: 



Anas ( ml \ , ■ 2 4vr^"s , 
-Cp— — 1 h... ]+lCF—mn^ (1 + 



1 1- . . . -|- i U-i? lllr)- 

k^ \ k^ J \k\ P 



= -Cp —asrmjj + ... 

r 2 

+i^asr^Tml + 7^ + Invr + ln(r /i)^ - 1^ + ... . (114) 

The dots stand for higher-order real and imaginary terms. In the Coulomb part, we have 
displayed only the leading term. In the imaginary part, the divergence comes from the 



Fourier transform of 1/|A;|^, which, in d dimensions, may be found in 3J|. Because we 
expand the gluon propagator in ml/k"^ the integral corresponding to the static quark self 
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energy has no scale and the matching gives 5m = 0. 

Next, we integrate out the scale mo- At one- loop level this corresponds to evaluating the 
contribution to the potential of the diagram generated by the singlet-octet vertex (dipole 
interaction) in the effective theory, which is shown in Fig. \5\ the HTL resummed gluon prop- 
agator is used. Therefore, the contribution is the same as the one calculated in Sec. IIVD 21 
and given in Eq. (!87|) . Summing its diagonal element with Eq. (11141) gives 

Vs{r) = -Cf— —asrmjj + —asr + . . . 

-i^a,r'Tml(^-2jE-HrmDf + l^+... • (115) 

We see that in the sum the divergences of Eqs. flll4p and (1871) cancel each other providing a 

Cf 

finite physical result. The term — as r'^m^ in the real part is suppressed by a factor rrriD 

6 

Cf 

with respect to — — ag r and will be neglected in the following. Note the appearance of 

the logarithm InirmD)'^ opposed to the appearance of the logarithm In — tt- in Eq. (1^21) . In 
the first case, the logarithm signals that divergences have been canceled when integrating 
out the scales 1/r and mo, in the second case it signals that divergences have been canceled 
when integrating out the scales T and m^i. The real and imaginary parts of Eq. (IllSp 
can be also obtained by expansion in rm^) of Eq. (11121) and — r/2, as defined in Eq. (11131) . 
respectively. 



F. Singlet static energy for 1/r » rriD 

The real part of Eq. (IllSp provides the static quark-antiquark energy for 1/r ^ mr,, 
whose leading thermal contribution is 

Cf 

5Es = — —asrm\,, (116) 

and minus twice the imaginary part of Eq. (IllSp provides the static quark-antiquark thermal 
decay width 



V = ^asr'Tm'j,\^-2-iE-HrmDY + -j ■ (117) 

Also in this case the singlet to octet break up process provides a contribution to the thermal 
width, which is 5T = -NICfoiIT. This contribution is parametrically suppressed by a 



factor {/W/mnY with respect to the one in (11170 . However, depending on and rm^j, it 
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may still contribute to a large fraction of F, see Fig. [161 Note that the condition mo ^ 
requires as(l/r) ^ IjNc x rniD to be fulfilled. 




as = 0.2 



0.2 0.4 0.6 0.8 1 

rniD 



FIG. 16: T /5T vs rniD for different values of a^. 



VII. CONCLUSIONS 



We have studied the real-time evolution of a static quark- ant iquark pair in a medium of 
gluons and light quarks characterized by a temperature T. We have addressed the problem 
of defining and deriving the potential between the two static sources, and of calculating 
their energy and thermal decay width. In the different ranges of temperature considered, we 
have set up and worked out a suitable sequence of effective field theories. Our framework 
has been very close to the modern EFT treatment of non-relativistic and static bound 
zero temperature [lC|], but complicated by the existence of the thermal scales T 
and tud- We have assumed that all the energy scales are perturbative and worked in a 
strict weak-coupling framework. This had two consequences: first, we could exploit the 
hierarchy T ^ second, the potential that we obtained is valid in the short range. In 
an EFT framework, the potential is the r-dependent matching coefficient that appears in 
front of the four-fermion operator that destroys and creates the bound state, after having 
integrated out all scales above the bound-state energy. Higher-order operators give lower 
energy contributions, entering into the computation of physical observables, but not in the 
Schrodinger equation that governs the motion of the bound state and hence are not of a 
potential type. 
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A systematic treatment of non-relativistic bound states in a thermal medium, in an EFT 
framework and in real-time formalism, has to our knowledge not been presented so far. We 
have devoted several parts of this work to set up a proper real-time formalism for static 
sources. The main outcome of this more formal aspect of the paper is in Eq. fIBUl) . which 
expresses the real-time quark-antiquark propagator as an infinite sum of free propagators 
and potential or mass-shift insertions. In all the considered dynamical regimes, the structure 
of the potential is such to satisfy this equation, see Eqs. fl58l) . fl9T|) . fllOOp and fllOip . 

We have considered a wide range of temperatures and provided the leading thermal 
effects to the potential. The results may be summarized in the following way. (1) If the 
temperature is smaller than or as large as the bound-state energy scale, there are no thermal 
contributions to the potential, which is simply the Coulomb potential fISS]) . Thermal effects 
enter in physical observables in loop corrections induced by low-energy gluons. The static 
energy and the thermal width of the system have been calculated in ( IHHj) and (!89|) . The 
system exhibits a thermal width, because, due to thermal fluctuations, at short range a color 
singlet quark-antiquark state may break up into an octet state and gluons. This is a novel 
feature of quark-antiquark bound states in a thermal medium that, to our knowledge, has 
not been considered quantitatively before. The singlet to octet break-up mechanism provides 
the dominant contribution to the thermal width when the temperature is as large as the 
bound-state energy. (2) If the temperature is larger than the bound-state energy scale but 
smaller than 1/r and if m^) is smaller than or as large as the bound-state energy scale, then 
the leading thermal correction to the Coulomb potential is given by the sum of Eqs. ( 1731) and 
( I82l) . The potential develops a real and an imaginary part. Both the singlet to octet thermal 
break-up mechanism and the imaginary part of the gluon self energy induced by the Landau 
damping phenomenon contribute to the imaginary part. If mo is smaller than the bound- 
state energy scale then the singlet to octet thermal break-up provides parametrically the 
dominant contribution to the decay width. The imaginary part of the potential originated 
by the imaginary part of the gluon self energy is divergent and needs to be regularized. In 
physical observables, the divergent part cancels against contributions coming from the scale 
mc). (3) If also is larger than the bound-state energy scale, then HTL of the type shown 
in Fig. [5] contribute to the potential as well. The expression of the potential turns out to 
be finite at the considered order and is given by Eq. fl92l) . Equation fl92l) is finite because it 
provides the leading thermal correction to the static energy and to the decay width in the 
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expansions in AV^/T, mo/T and rm^- The same result has been also derived in perturbative 
QCD by calculating the Fourier transform of the diagrams in Fig. [9] and Fig. [TTJ (4) If the 
temperature is larger than 1/r but m£, is smaller than 1/r, the static potential is given by 
Eq. (11151) . If is also smaller than or of the same order as AV then the potential is 
given by Eq. flll4p and divergences cancel in physical observables against loop corrections 
from the scale rriD. (5) Finally, if m^i is of the order of 1/r the static potential is given 
by Eq. (11 111) : this result agrees with the earlier finding of 5|. Equations (11151) and (II lip 
are finite because, in the kinematical regions of validity, they provide the leading thermal 
correction to the static energy and the decay width (see Eqs. (I116I) .( |TT71) . (11121) and (11131) ). 
In the temperature ranges (3), (4) and (5), the parametrically dominant contribution to 
the thermal width comes from the imaginary part of the gluon self energy. 

While this work was being carried out and completed some papers appeared dealing 
with some of the issues addressed here. In [35|, l36|, among others, the role of the gluon 



condensate was considered. The gluon condensate enters the expression of the mass of the 
quarkonium if the typical binding energy is much larger than the hadronic scale Aqcd- The 
potential is therefore not affected, while the effect of the gluon condensate on the mass is 
parametrically smaller than the leading relativistic corrections. If the temperature scale is 
larger than the binding energy then the temperature enters the potential and the effects 
have been described in the present work. If T ^ 1/r, temperature effects are carried by the 
dipole-dipole interaction shown in Fig. [5] and corrections to it. Only when the time scale 
of the chromoelectric correlator is smaller than the binding-energy scale, the chromoelectric 
correlator reduces to the local condensate (see Sec. IIVCI) . Low temperatures, smaller than 
the binding energy, affect the condensate but not the potential, which remains Coulombic. 
In this case, thermal effects (in the static limit, they have been evaluated in Eq. fl90|) ) are 
parametrically smaller than the leading relativistic corrections. It remains to be clarified if 
a temperature scale below the binding energy is above the critical temperature: this may 
depend on the fact that the temperature scale is T or a multiple of ttT. In js^], some of 
the issues discussed here about the potential in real-time formalism have been addressed. 
In particular, static particles in real-time formalism have been considered and Eq. (Illip has 
been derived for a hot QED plasma in a real-time framework. It has been also pointed out 
that the real part of the static potential (II IIP agrees with the singlet free energy and that 
for large separations r ^ l/rriD, T/2, as defined in Eq. (I113p . gives twice the damping rate 
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of an infinitely heavy quark. Finally, in 38|] a study of non-relativistic bound states in a hot 
QED plasma has been carried out in a non-relativistic EFT framework, which is similar to 
the one presented here. 

There are many possible developments of this work some of which have already been 
mentioned in the previous pages. Here we stress few of them. First, the construction of a 
full EFT for non-relativistic bound states at finite temperature requires to be completed. 
We have focused in the paper on the quark-antiquark color singlet state, but a complete 
identification and study of all relevant degrees of freedom that appear once the thermal 
energy scales have been integrated out. ,s stU. to be done. Th.s ,nay reqn,,. the constrnetion 
of an EFT that includes the dynamics of gauge fields below the scale m^j [3^ . Second, in the 
EFT framework presented here, the study of quark-antiquark states at large but finite mass, 
i.e. actual quarkonium in a thermal medium, should be addressed. As argued along the 
paper, the static limit provides the first piece of a 1/m expansion; higher-order corrections 
may be systematically implemented in the framework of NRQCD and pNRQCD. Finally, 
although the short- distance analysis performed in this work may provide a valuable tool 
for studying the thermal dissociation of the lowest quarkonium resonances, the inclusion in 
the analysis of the non-perturbative scale Aqcd inay become necessary for studying excited 
states. 
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